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Abstract 


This paper showed that Poisson brackets in quaternion variables 
can be obtained directly from canonical Poisson brackets on cotan¬ 
gent bundle of SE( 3) (or 50(3)) endowed by canonical symplectic 
geometry. 

Quaternion parameters in our case are just dynamic variables in 
canonical Hamiltonian mechanics of a rigid body on T*SE{ 3) 

The obtained results based on quaternions representation as ex¬ 
plicit functions of rotation matrix elements of 50(3) group. 

The relation of obtained Poisson structure to the canonical Poisson 
and symplectic structures on T*S 3 were investigated. 

To derive the motion equations of Hamiltonian dynamics in quater- 
nionic variables it is proposed to use the mixed frame of reference 
where translational degrees of freedom describes in the inertial frame 
of reference and degree of rotational freedom in the body frame. 

It turns out that motion equations for system with Hamiltonian 
of the rather general form can be written in algebraic operations on 
quaternions. 

Keywords: quaternion, symplectic structure, Poisson structure, 
Lie-Poisson brackets, Liouville one-form, SE(3), S0(3). 
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1 Introduction 


This article is describing main results of Hamiltonian dynamics of a rigid 
body in quaternion variables. 

There are many of publications that are devoted to application of quater¬ 
nions in mechanics of a rigid body. The most of them refer to kinematics of 
a rigid body, i.e. description of rigid body orientation in quaternion param¬ 
eters [1]. 

Hereinafter we specify those rare articles that were devoted to description 
of dynamics with using quaternions. Kozlov’s work [2] gives rigid body dy¬ 
namics in Lagrangian representation based on classic approach developed by 
Poincare. Borisov and Mamaev [3, 4, 5] proceeding from the deep relations 
between quaternion algebra and groups SO( 3), -SO(4) have proposed Poisson 
brackets between parameters of quaternion and components of intrinsic an¬ 
gular momentum of a rigid body and considering them as some generators of 
a Lie-Poisson structure. Therewith condition of normalization of quaternion 
per unit that is required for a rigid body description emerges as the special 
value of the Casimir function in that Lie-Poisson structure. 

Expression of matrix element of the group in quaternion parameters is well 
known [1, 6, 4, 5] that is also clearly demonstrates structure of orthogonal 
matrices. Per se inverse problem of expressing quaternion parameters in 
terms of elements of the corresponding rotation matrix is not so difficult. 
But in literature we have found (post factum) only one article devoted to 
solution of that task [7]. Main goal of this short but contensive article is 
the comparative analysis of necessary number of operations for computation 
through matrix or quaternion. 

In terms of quantity was proven that quaternions computation is more 
effective. Using of quaternion also has advantage of stability for numerical 
integration of motion equations (see [4, p.104]). 

From general relations given in work [8] one can deduce following expres¬ 
sions for Poisson brackets between elements of rotation matrix and compo¬ 
nents of intrinsic angular momentum in inertial frame of reference. 

Here and further for determinacy will consider SE( 3) that describes not 
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only rotational but also translational degrees of freedom of a rigid body. 


( 1 ) 


'{xi,Xj} = 0 ; 

{ x iiPj} — fiiji 

{ x ii Qjk} 0 , 

^{xi,mj} = 0 ; 


{PiiPj} 0) Qjk} E-ijlQlk) 
{Pi, Qjk} = 0; {mi, rrij } = £ijim h 

{Pi, mj} = 0 ; 

}Qij i Qkl} o, 


where Xi - coordinates of the body center of mass, pi - momentum compo¬ 
nents of translational motion, Qjk - elements of rotation matrix that describe 
the body orientation with respect to inertial frame of reference, rrij - compo¬ 
nents of intrinsic angular moment of the body with respect to inertial frame 
of reference. 

In fact expression ( 1 ) will also be proved during the explanation by the 
different ways. 

From the explicit expression of quaternion parameters as functions of 
elements of rotation matrix from (1) one can deduce required Poisson brackets 
between components of quaternion and angular momentum components of 
the body. Obtaining of these relations is main goal of the paper. 

It turns out that expressions obtained in this way have the same form 
that previously discussed Lie-Poisson brackets between generators of Poisson 
structure mentioned above (see [4, (2.7), p. 103]). 

But obtained relations have other core idea. Quaternion parameters in 
our case are not generators of Poisson structure, but dynamic variables in 
canonical Hamiltonian mechanics on T*SE( 3) (it is not Lie-Poisson struc¬ 
ture). In our case structural tensor of Poisson brackets of Lie algebra is 
not nondegenerate, and, so, Casimir functions are absent and normalization 
condition for quaternion just reflect relations between quaternion dynamic 
variables on group. To illustrate that fact let us consider example. Let is 
one of Euler angel. Then dynamic variable cos 2 ( 99 ) + sin 2 ( 99 ) is identically 
equal to one, but it is not the Casimir function. 

After obtaining Poisson structure for quaternion variables as dynamic 
variables for a rigid body (on T*SO( 3)) we investigate its relation with 
canonical Poisson and symplectic structures on T*S 3 (the relation is local 
isomorphism). 

To derive equations of motion in Hamiltonian dynamics in quaternion 
variables is proposed to use a mixed frame of reference where translational 
degrees of freedom describe in inertial frame of reference, and the rotational 
in the body frame (a reference frame associated with a body). It turns 
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out that motion equations can be written through algebraic operations in 
quaternions. 

Symplectic and Poisson structures on T*S 3 can be set completely inde¬ 
pendently from algebraic structure of quaternion variables, but one can see 
that Poisson brackets for these dynamic variables is completely determine 
laws of multiplication in quaternion algebra. 

Hamiltonian dynamics of a rigid body in quaternion variables is an in¬ 
teresting point of intersection of two pioneer directions in mathematics that 
were offered by Hamilton: powerful stream of Hamiltonian formalism that 
covering almost all areas of mathematical physics and not so strong but sur¬ 
prisingly beautiful stream of quaternionic applications in various fields of 
mathematics and physics. 
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2 


Hamiltonian formalism on T*SO(3) 


2.1 Representation of right trivialization of T*SO{ 3) 

Let’s bring are some useful relations for group SO(3) and its cotangent bundle 
many of which can be found in [9, 10]. 

Group 50(3) is formed from orthogonal unimodular matrices R i.e. R 7 = 
R 1 , det(R) = 1. Accordingly the Lie algebra so( 3) is formed from anti¬ 
symmetric 3 x 3-matrices with a Lie bracket in form of matrix commutator. 
Let’s introduce isomorphism of vector spaces": M 3 — > so( 3) such that [10, p. 
285]: £ kl = —euati, & = -\eikHw where e m — Levi-Civita symbol. Then 

!?*? = £ X 77; [£, r)\ = £rj - f/? = £ x rj; 

\(^r7) = -itr(^); B^B 1 = B| 

The scalar product above allows to set equivalence of a Lie algebra and its 
dual space so(3) ~ so(3). Symbol mean a diffeomorphism. As a rule, 
diffeomorphisms used below have a simple group-theoretical or differential- 
geometrical meaning that is explained in cited literature. 

The representation of right trivialization corresponds to inertial frame of 
reference that is shown in [9, p. 314], we have 
TSO{ 3) ~ 50(3) x so(3), T*SO{3) ~ 50(3) x so( 3)*. 

Then right and left group action of 50(3) on T*SO( 3) (Cotangent Lift 
[10, p. 166]) takes the form 


Ag : (R, 7r) G T*SO( 3) ->• (RB : tt), Be 50(3); 

Lg : (R, tt) G T*50(3) ->• (BR. BlrB 3 ) = (BR,Adg-i7r) 


2.2 Symplectic and Poisson structures on T*SO( 3) 

The Liouville 1-form on T*SO( 3) ~ 50(3) x so(3)* is given by: 

= -|tr(»Wft) = (3) 

where = (dR)R _1 — right-invariant of the Maurer-Cartan 1-form. 

Then for canonical symplectic 2-form, using Maurer-Cartan equation [10, 
. 276] we have: 

nZa SO(:i) = -de T * so{3) = 5R 1 a d-Ki - m[SR, dRy = ( 4 ) 

1 1 

^£ijkdRjk A djlj T ~Tli£ijkdRjs A SR s k 

Any given symplectic structure defines a Poisson structure on the same 
manifold: 

{F,G}(z) = n(fr(z),fc(*)) = d^F = -d^G (5) 

where for vector field the relation i^ a G, = dG is fulfilled. 

Let’s consider the elements of the matrix R and components of the mo¬ 
ment 7T as the dynamic variables on T*SO( 3) and we obtain the following set 
of Poisson brackets that completely determine Poisson structure on T*SO( 3): 

{RijiRkl} 0, {4Tj, -Rjfcj - EijiRtk, {7Tj, 7Tj} • (6) 

Notice that in inertial frame of reference Poisson brackets for matrix 
elements R is grouped together in columns, for example, Poisson brackets of 
3-rd column expressed only through elements of 3-rd column. 
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3 Quaternion algebra 

Quaternions form an associative algebra with unit eo and generators 
e*, i = 1, 2,3 that satisfy to defining relation 

C r C s 0 ^rst&t ( 1 ) 


or, that is equivalent to 


[e r ,e s ] e r e s e s e r C 2s rs ^ci^ 
— 2(5t*sCo* 


( 2 ) 


Thus quaternions form a 4-dimensional vector space over the field of real 
numbers 

q = q°e 0 + q 1 e l + q 2 e 2 + q 3 e 3 (3) 

or, through components of a 4-dimensional column vector: q = (g°, q 1 , q 2 , q 3 ). 

The component q° is a scalar part of a quaternion q and components 
q 3 ,q 2 ,q 3 grouped as a vector part. Thus quaternion can be written as q — 

(q°,q)- 

If q° = 0 then permissible form q = q, and such quaternion called pure 
quaternion [6, p. 301]. The pure quaternions form a linear subspace of 
quaternions algebra, but not subalgebra, because associative product of two 
pure quaternions leads to quaternion of general form 

From (1) follows law of quaternion multiplication a and b 

ab = (a°e 0 + a r e r )(b°e 0 + b s e s ) = (a°e 0 + a)(6°e 0 + b) (4) 

= (a°6° — (a, b))e 0 + a°b + b°a + a x b 

Then for pure quaternions through law of associative multiplication we 
have the following expressions of scalar and vector product. 


I (x,y) = —(xy + yx)-, 

\x x y = \(xy - yx) 

Operation of quaternions conjugation is introduced as follows 

eo f e,t 

= {q°, -q) 


(5) 

( 6 ) 
(7) 
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Thus pure quaternions is completely characterized by the following prop¬ 
erty 

x1 = —x (8) 

From (7) follows 

qq< = q'q = (q°q° + (q, q))e 0 = (q°) 2 + (q 1 ) 2 + (q 2 ) 2 + (q 3 ) 2 (9) 


Thus, from (4,7) follows 

(aby = b^a) 

Let introduce norm of quaternion 

kl = V (qq 1 ) 


( 10 ) 

( 11 ) 


then 

\ab\ = |a||6| 


( 12 ) 


In addition from (9,11) follow a simple representation of inverse quater¬ 
nion 


-i _ 

M 2 


(13) 


Formulas (9) and (13) indicate that all quaternions but excluding a zeroth 
have its own inverse, i.e. algebra of quaternions is the division ring. 

From (12) and (13) that a unit quaternions, i.e. quaternions that norms 
are equal to unit form a group. 
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4 Presentation of rotations in quaternions 

Let’s consider linear subspace of pure quaternions as 3-dimension Euclidean 
space with inner product that determined by formula (5) §3. 

Let’s show that inner automorphism of quaternion algebra that generated 
by unit quaternion q, translate space of pure quaternions into itself. 

The law of transformation can be written as 

(1) x = qxq~ l = qxq^ = Q[x\, \q\ = 1. 

Let’s show that resulting quaternion x also is a pure quaternion, indeed, 


(qxq^)^ = (g i ')Vg' 1 ' = —qxqK 


( 2 ) 


so, operator Q[x\ is a linear operator that acting in a subspace of pure quater¬ 
nions. 

Let’s show that operator Q[x] preserves scalar product of vectors, using 
expression of scalar product through associative multiplication from (5) §3 



Thus, following relations for invariance of operator Q are correct (proved 
similarly) 



(Q[x\,Q[y)) = (x, y)] 

Q[x] x Q[y\ = Q[x x y}; 
(Q[z\,Q[x\ x Q[y}) = (z,x x y) 


(3) 


The first relation (3) means that operator Q is orthogonal and third 
relation that it is unimodular, i.e. Q is proper rotation (Q 6 SO(3)). 

The particular interest in context of posed a problem is explicit form 
of matrix elements of operator Q that can be obtained from (1) by using 
multiplication law (4) §3 
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that for unit quaternion is equivalent to following expression for matrix of 
operator 


(5) Q 


4 + ?i ~ 4 ~ Qi 

2(gi<?2 + 90*73) 
2(<M3 - qoQ2) 


2(9192 - <?o<?3) 

qi - qi + 4 - qi 

2(9293 + 9o9l) 


2(9193 + 9092) 
2(9293 - 9o9 1) 

4 -qi- 4 + 4. 


Let’s show that q —> Q is mapping on whole SO(3) group. 

In particular case when qi = 0, g 2 = 0 this matrix takes the form 


Q 


4 - qi -2<ms o 
2<ms qi - qi o 

0 0 ql + ql 


If q 0 = cos q 3 = sin then 


Q 


cos 6 — sin 6 0 
sin 6 cos 6 0 

0 0 1 


This is a matrix of rotation around Z axis through an angle 6 (counterclock¬ 
wise). 

Likewise, we will get rotations around x and y axes. 

Notice. It is useful to note that an arbitrary unit quaternion can be 
represented as 


( 6 ) 


q = cos [ipj 2)e 0 + sin (<^/2)e, 


where e is pure and unit quaternion that sets the rotation axis and (f is 
corresponding angle of rotation around that axis. Expression (6) that is 
considered as a function ip is one-parameter subgroup of the group of unit 
quaternions. 

Because of T : q —» Q is a homomorphism, so product of such rota¬ 
tions be contained in the image of this homomorphism. It is known that any 
rotation can be obtained as a product of rotations around the axes of Carte¬ 
sian coordinate system. So, matrices (4) and (5) are record of an arbitrary 
element of the group SO(3) in terms of quaternion parameters. 

Thus formula (1) defines a homomorphism of the group of unit quater¬ 
nions on SO(3). At the same time, these groups are locally isomorphic and 
the group of unit quaternions is double covering group SO (3) [6]. Indeed, 
from (1) it is obvious that quaternion q and (—q) gives the same rotation Q. 
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5 Quaternion representation in terms of ro¬ 
tations 

In previous section it was shown that every rotation SO (3) corresponds to 
2 and only 2 unit quaternions of the different sign. There is the problem of 
explicit, expression of this functional dependence. 

Matrix Q that corresponds to quaternion q has following elements 

(1) Qik (2go T 2 q%qk ^QoQj^ jik 

Let’s obtain the trace of matrix with taking into account that q is a unit 
quaternion 


Sp(Q) = 4qjj - 1 


( 2 ) 


i.e. 


Qo = J(Sp(Q) + 1) 


( 3 ) 


The antisymmetric part of matrix Q has a simple form, so 

QoQi ~^ijkQ jk 


( 4 ) 


Thus where go 7^ 0 one can explicitly express components of quaternion 
corresponding to a given rotation matrix, through its elements 




There are two solutions to this system, corresponds to two various choices 
of sign of root in the expression for go. 

If we want to have a solution of these functional equations in the neighbor¬ 
hood of quaternion with zero scalar part (Sp(Q) = —1), i.e. in neighbourhood 
of pure quaternion then formulas (5) are not suitable. 

It should be also mentioned that for numerical computations the difficul¬ 
ties can also arise for small but non-zero values of g°. Therefore E. Salamin 
who investigated the problem of comparative evaluation of effectiveness of 
numerical computations with using orthogonal matrices and quaternions [7], 
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suggested to use next set of formulas to find components of quaternion that 
follows from matrix Q 


( 6 ) 


P) 


Qo = i(l + Q n + Qn + Q33); 
Qi — |(1 + Qu ~ Q22 ~ Q33)', 
( li = |(1 — Qu + Q22 ~ Q33); 


QoQi 

qoQ2 

qoq3 

9 i ?3 

l q-iqi 


Qn — 

Q 

22 + C 

\(Q?,2 

— 

Q23)', 

i (Ql 3 

— 

Q31); 

KQ21 

— 

Q12)', 

\{Ql2 

+ 

Q21)', 

\(Ql 3 

+ 

Q31)', 

i (Q23 

+ 

Q32 ); 


For numerical computations it is optimal to choose of the maximum com¬ 
ponent of quaternion (by its absolute value) from (6) (choice of the sign of 
component will determine the sign of quaternion), and then the others com¬ 
ponents can be found from (7). 

It is clear that for unit quaternion at least one of the component will be 
different from 0. 

Note that solution (5) that found above corresponds to 1-st line of (6) 
and first 3 lines in (7). 
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6 Poisson brackets with components of quater¬ 
nion 

As it is follows from the foregoing description quaternion variables can be 
regarded as dynamic variables in the phase space of a rigid body. The only 
difference from usual dynamic variables in two-valuedness of quaternion vari¬ 
ables. 

This is not an obstacle to use of quaternion variables to describe of dy¬ 
namics of a rigid body. By selecting a certain sign of these variables at the 
initial point of the trajectory in phase space we choose the sign in accord 
with the continuity throughout the trajectory. 

Also, if we restrict some neighborhood of a point of phase space of a rigid 
body then it is possible to choose one of branch of expressions in quaternion 
variables express in terms of elements of rotation matrix. 

To derive Poisson brackets with components of quaternions we need to 


com P ute Si and S 


For instance, 



where Poisson bracket {7T ? ;, Q k i} is known from (1) §1. 
From (5) of previous section we have 




If we take into account { 71 Q ki } = £i kn Q n i with (2) we get 


( 3 ) 



and 


( 4 ) 




Since {Qij, Q k i} = 0 we finally obtain 

3 ; 
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As mentioned in the introduction, such Poisson brackets were obtained 
by Borisov and Mamaev as an expressions for generators of some Poisson 
structure of Lie-Poisson. In addition the value of C(q) = ( q°) 2 + ( q 1 ) 2 + 
( q 2 ) 2 + ( q 3 ) 2 is Casimir function in Poisson structure and so in order to get 
description of dynamics of a rigid body we should pass on to a symplectic 
leaf of this Poisson structure that corresponds to the value C(q) = 1. 

6.1 The properties of invariance of Poisson structure 

on T*(SE( 3)) (T*(SO( 3))) 

Take note on Poisson brackets (1) §1 that contain matrix elements of rotation. 
For example, look at expression of 

(1) Qjkf E-ijlQlk 

lets multiply it on fixed matrix B € SO (3) from right side. Since elements 
of the matrix B are constants then we can put them into Poisson brackets 
in the left side of (1) then the expression is valid 

(2) {v( QB ) jn } = £iji(QB)i n 

If we performs the same conversion with all such Poisson brackets we 
can see that matrix elements of rotation P = QB will satisfies to all those 
Poisson brackets as matrix elements of initial rotation Q. 

Let’s show that Poisson brackets of previous section can be written in a 
more compact form by using law of quaternion multiplication. 

Let’s multiply quaternion q = q 0 e o + q^k from left side on e* 

Cj*? qo&i T qk(&i&k) Qo&i T qk( 0 T ftC0 "L T qk^ikj^-j 


i.e. 

( 3 ) etq = —qieo + q^S^ej + qk£ikj£j 


or 


( 4 ) 


(e*g) o = -g»; 

(AA)i &ijkqk T q(A ]: 
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(5) 


Then 


{<?o,7Ti} = | (e^o; 

{/U, t} = |(eig)j-; 


i.e. 


( 6 ) 


{%.,"■(} = A = 0,1,2,3; 


Formally one can multiply (6) from left side on fixed quaternions eo, ei, e 2 , e 3 
and then summing by //-index. In addition using as constant we can insert 
them into Poisson bracket. 

The result is 


(7) 


{7Ti,q} = ~-e t q ; 


In the same way as for matrices when we multiply (7) from the right 
side on a some fixed quaternion b,\b\ = 1 we can see that new quaternion 
dynamic variable p = qb satisfies same Poisson brackets (7) like q. In fact 
it is a simple algebraic transformation is equivalent to applying matrices of 
right translation (2) §2 to the set of relationships (7). 

Therefore p satisfies to all Poisson brackets that are analogous to Poisson 
brackets in (5) of previous section. 

6.2 Poisson brackets with quaternion components, the 
special case 

Strictly speaking the obtaining relations (5) §6 that are present above will not 
be correct in the neighbourhood of pure quaternion Sp(Q) = —1 —> q 0 = 0, 
because (5) §5 is not feasible the solution at this point. 

But this not means that formulas (5) §6 in the neighbourhood of pure 
quaternion is not valid. A quaternion with go = 0 does not looks like a 
critical point for these formulas. However the proof of these formulas in 
general requires other approaches. 

First of all one can use Salamins solutions of (6,7) §5 by taking as a 
pivotal the component of quaternion that is not zero in this neighborhood, for 
example q 1 . Unfortunately such direct approach leads to very cumbersome 
and non-transparent computations. 
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Let’s consider properties of invariance that were listed in previous section, 
as well as the fact that the map q —y Q is a homomorphism on group SO (3) 
as it was described above. 

Let’s suppose that we need to obtain Poisson bracket with quaternions 
in the neighbourhood of pure quaternion q with go = 0. 

Lets consider such fixed quaternion b : \b\ = 1 that dynamic quaternion 
variable p = qb has value p° ^ 0 at this point (this always can be done for 
q : \q\ = 1 as it was shown above). Quaternion p mapping into the matrix 
P G SO (3) for the chosen homomorphism (p —y P ). 

Quaternion variable p has same relation to matrix-variable P just as the 
relation of variable q to variable Q when deriving Poisson brackets with 
quaternions that was described above and moreover p° ^ 0. 

Consequently for quaternion variable p Poisson brackets (5) §6 are ful¬ 
filled. By applying the inverse transform q = pb~ l to the obtained formulas 
we will got (5) §6 but for quaternion variable q. 
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7 


Tangent and cotangent bundle over S 3 


7.1 Tangent bundle to group S 3 C H 

Since group S 3 is a subset of associative algebra HI (so and linear space) then 
a curve on the group is also a curve in the linear space. Therefore the vector 
of speed £ of this curve i.e. element £ € T(S 3 ) we can considered as the 
element of HI. 

(1) q(t)q(ty = 1— >qq ] +qq ] = (q,q) = 0 

So, tangent vectors to S 3 in the point q can be presented as quaternions 
£ that are orthogonal to quaternion q. 

The left (right) translation of tangent vector in this presentation coin¬ 
cides with product of quaternion that are tangent vector on quaternion that 
generates left (right) translation. 

Really, let q(t)\ t =o = a 

(2) T a L b ■ q = ^(L b q(t))\t=o = ^(bq(t))\ t=0 = bq 

Proposition 1. If e — pure quaternion then £ = L q e = qe e T(S' 3 ) 
and, vice versa, if £ 6 T(S 3 ) then g -1 £ = gt£ — pure quaternion (q 6 S 3 ). 
Similar statements are correct and for right translations. 

□ 

The fact that e - pure quaternion can be express in terms of 
e + e 1 " = 0 —» eel + e o et = 0 —> (eo, e) = 0 —> {L q e 0 , L q e) = 0 
and, vice versa 

(Q, 0 = 0 —» (Lq-iq, L q - 1 £) = 0 —> (e 0 , g _1 £) = 0 —> g^ x £ + (g _1 £) f = 0 

■ 

Conclusion. The tangent space at the group identity can be identified 
with subspace of pure quaternions. 

On T(S 3 ) from HI the Euclidean metric is induced. In particular T e (S 3 ) 
is a subspace of pure quaternions and at the same time it is 3-dimensional 
Euclidean space. 
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7.2 The Lie algebra of S 3 C IHI 


As shown above an arbitrary unit quaternion can be represented as cos (ip/ 2)eo+ 
sin(<p/2)e, where e - pure and unit quaternion. 

Chosen the one-half angle as a parameter of the one-parameter subgroup 
connected with the interpretation of homomorphism T : S 3 K > SO (3) wherein 
e defines an axis of rotation and ip - corresponding the angle of rotation 
around this axis. 

For study of the Lie algebra Lie(S 3 ) will be more convenient to present 
one-parameter subgroup in form 

(1) exp (te) = cos(t)eo + sin(t)e 

Accordingly the pure unit quaternions e can be considered as an element 
of Lie (S' 3 ) algebra. 

It is reasonable to assume as a basis for Lie (S' 3 ) set of generators e*, 


i = 1,2,3. 

Then formula (2) §3 provides structure constants of this Lie algebra 



( 2 ) 


So, from the above it follows that Lie(S' 3 ) = T e (S 3 ) can be identified with 
subspace of pure quaternions. 

Proposition 2. Lie(S' 3 ) = T e (S 3 ) subspace is revealed an Euclidean 
space where operator Ad q G S'0(3). 

□ 

Indeed, 


(3) Ad q [e] = — (gexp(ie)g = qeq 


,-i 


qeq ] = T(q)[e] 


Now it is easy to derive expression for Lie bracket in quaternions 



(3a) 


ad^r/] =£r]-ri€ = [£, rj] = 2£ x ry 
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7.3 The dual space to Lie algebra and bundle T*(5 3 ) 

The metric on Lie(S 3 ) = T e (S 3 ) allows us to identify Lie*(S 3 ) and Lie{S 3 ). 

Indeed, pure quaternion n defines linear form on Lie(S 3 ) = T e (S 3 ) by the 
formula 

(1) kK] = <k, 0, v?eL S e(S 3 ) 

Let’s compute the operator Ad*[/x] 

(2) (Ad*[/z], £) = (m, Ad,[£]) = (m, q^q~ x ) = {q~^q, 0 = (gW £) 

or 

(2a) Ad*[/j] = q~ 3 nq = q ] fiq 

Operator adjj can be found either as derivative with respect to t for Ad* [fj] 
where q{t) = exp(t£) or as operator adjoint to ad^. 

In any case we get the following result 

(3) ad£[/*] = [/*,£] =2/z x £ 

In the light of above said and bearing in mind the existence of 2-directional 
invariant metrics on S 3 we can identify T*(S 3 ) and T(S 3 ). 

7.4 Algebra isomorphism Lie(S 3 ) and Lie{SO{ 3)) 

Let’s consider expression of right-invariant Maurer-Cartan 1-form that was 
introduced in subsection 2.2. and represent it through quaternion variables. 

Proposition 3. The following formulas are valid. 


(1) {dQQ 1 ) ifc = -2 ( Qidq k - q k dqi + ( q 0 dq r - q r dq 0 ) e ikr ) 


( 2 ) 


(dQQ )j y^irsidQQ )rs 

= 2 ( q 0 dqi - q { dq {D + £ irs q r dq s ) = 2(dqq~ 1 ) i 
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□ 

We have 

(3) Qij (2^0 l)^jj 4" ^QiQj ^QoQr^-ijr 

(4) dQij = 4q 0 dq 0 5ij + 2q l dq 3 + 2q J dq i - 2q 0 dq r £ ijr - 2q r dq 0 £ ijr 

(5) Qj k = Qj k = ( 2( ?o - + 2 q 3 q k + 2q 0 q s e jks 

From here by using ql + q 2 = 1 —> qodqo + (q, dq) = 0 

( dC)C) j jk (~^qi)dqi)dij ~\~ 2qidqj *h 2qjdqi 2qodq r £ij r 2q r dqaEjj,- j 
x ((2 ql - 1 )S jk + 2q 3 q k + 2q 0 q s £ jks ) 

= 4:(2q^-l)q 0 dq 0 S ik +2(2q^-l)q i dq k +2(2q^-l)q k dq i -2(2ql-l)q 0 dq r £ ik r-2(2ql-l)q r dqo£ikr 
+8q 0 dq 0 qiq k + 4 qidq 3 q 3 q k + 4<? 2 q k dq t - Aq Q q 3 q k dq r £ ijr - Aq 3 q k q r dq 0 £ ijr 
+8q 0 q s qodq 0 £ iks + 4 q 0 qiq s £ jks dq 3 - kqlq s £.j ks £ijrdq r - Aq r q s £ jks £ ijr q Q dq Q 

= A(2q^-l)q 0 dq 0 S ik +2(2q^-l)q i dq k +2(2q^~l)q k dq i -2(2ql-l)q 0 dq r £ ikr -2(2ql-l)q r dqo£ikr 
+8q 0 dq 0 qi.q k + 4 qiq k q 3 dq 3 + 4q 2 q k dqi - 4 q 0 q 3 q k dq r £ ijr - 4:q 3 q k q r dq 0 £ i3r 
+8qlq s dq 0 £ iks + 4 q 0 qiq s £ jks dq 3 
rs $kr Sis)dq r + 4 q r q s (6 ik S rs - 5 kr 6 is )q 0 dq 0 

= 4(2go - l)q 0 dq 0 5 ik + 2(2 ql - 1 )qtdq k + 2(2 ql - 1 )q k dqi 

2(2g 0 X^qodq-rEikf 2(2^ \')q r dqo£i kr 

+8qiq k q 0 dq 0 + 4 qiq k q r dq r + 4 q 2 q k dqi - 4 q 0 qjq k dq r £ ijr - 4 q 3 q k q r dq 0 £ ijr 
+8qlq s dq 0 £ iks + 4 q 0 qiq s £ jks dq 3 
+4 (qlq r dq r + q 2 q 0 dq 0 ) - 4q k qiq 0 dq 0 ~ ^qlqdq k 

= 4 ((2<?o - l)go^o + [qlqrdq r + q 2 qodqo)) $ik 
+2 (( 2 ?o - 1 ) - 2 ^q) + 2 (( 2 g„ - 1 ) + 2 q 2 ) q k dqi 

2(2g 0 1 ^)qodq r £i kr 2(2 q^ \^jq r dqo£i kr 

+4q i q k q 0 dq () + 4 qiq k q r dq r - 4 q 0 q 3 q k dq r £ ijr 
+8qlq s dq 0 £ iks + 4 q 0 qiq s £ jks dq 3 
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= 4 (<7g (qodqo + q r dq r ) + (q% + q 2 )q 0 dq 0 - q 0 dq 0 ) 5 ik 
+2 ((2<?o - 1) - 2go) ftrig* + 2 ((2^ - 1) + 2q 2 ) q k dq t 
2(2g 0 t)qodq r £ikr 2 (2<^q \^q r dqo£i kr 

+4qiq k q 0 dq 0 + 4 qiq k q r dq r - 4q 0 qjq k dq r £ ijr 
+8qlq s dq 0 £ iks + 4 q 0 qiq s £j ks dqj 
= -2q t dq k + 2q k dq t 

4 q{) qii d(fr£/ kr 4 “ 2qodq r £i kr 2q r dqo£i kr 

+4 qlq r dq 0 £ ikr + 4q 0 q i q r £ jkr dq j - 4q 0 qjq k dq r £ij r 
= - 2 q t dq k + 2q k dq { - 4$(q 0 dq r - q r dq 0 )£ ikr 
42 qodq r £j ikr 4 " 2q r (lq()£j kr 4 " 4q$qj {q%£ k j r qk^ijr)dq r 

Let’s use the identity 


qi£kjr qk£ijr ~ 2 £iks£smn{ - qm£njr Qn£mjr) 


^E'iks^Qm^njr^ smn Qn^mjr^smn) 


2 £ iks\Qm{^js^mr $rs$jm) Qn\^js^nmr ^rs^jn) 
^iks^Qr&js Qj& rs ) 


Then 


(■dQQ 1 )« = - 2 q t dq k + 2 q k dqi - 4q 2 (q 0 dq r - q r dq 0 )£ ikr 

-\~2qodq r £i kr 4 - 2q r dqo£i kr 4 ~ 4qoqj£j /ks (^q r 5j S qj5 rs ')dq r 
= -2qidq k 4 - 2 q k dqi - 4ql{q 0 dq r - q r dq 0 )£i kr 
+2q 0 dq r £ ikr + 2q r dq 0 £ ikr + 4q 0 £ iks (q s q r dq r - q 2 dq s ) 

= -2 q t dq k + 2q k dqt 

-4q 0 (q 2 dq r + q 2 dq r - q r q 0 dq 0 + q r q 0 dq 0 )£ ikr 
-\~2qodq r £j /kr 4 ~ 2q r dqo£i kr 
= -2qidq k + 2q k dq { - 2 q 0 dq r £ ikr + 2 q r dq 0 £ ikr 
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Thus, it is proved the relation (1), i.e. 

(■ dQQ = -2 ( Qidq k - q k dqi + ( q 0 dq r - q r dq Q ) e ikr ) 

Going to the vector representation we obtain 

(dQQ )j ~£jik(dQQ )jfc e ji k (qidq k q k dqi^j T ( qod(/r qrdqo) £jik£ikr 

= 2 £ jik qidq k + 2 (g 0 dg r - g r dg 0 ) &jr 
= 2 (godgj - qjdq 0 ) + 2 £ jik qidq k 
Hence was proved the relation (2), i.e. 

(dQQ )j — £i rs (dQQ ) rs 2 (q'o dq% q%dqQ T £%rsqrdqs) 2(dqq )j 


Let’s call arithmetic vectors through bold symbols, i.e. sets of numbers 
that are vector components, and let for angular velocity and intrinsic angular 
momentum they should be components in inertial frame of reference, but for 
quaternions they should be in basis of generators of quaternion algebra. 

Then in these coordinates 


(6) dT e = 2 • Id 

(7) u> = 2dT c [£], ueLie(SO( 3)), £ e Lie(S 3 ) 

We have 

(>«[*] = !«: 

(8) < rfT.KI = iw'; 

LU[2£ x(l=ux a/; 
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8 Poisson structure for dynamic quaternion 
variables 

8.1 Jacobi identity for structure tensor 

Poisson structure for dynamic quaternion variables can be entered axiomat- 

ically. 


'Un,qv} = 0 , fi,v = 0 , 1 , 2 , 3 ; 

{jJ-ii \ Qi i 

{ l^'i • Qj\ E-ijkQk Qo^iji 

„ {/L; /b l 2 Eijifli, 


( 1 ) 


But for that we must prove the Jacobi identity for structure tensor of 
quaternionic Poisson brackets. 


Proposition 4. The structure tensor of quaternionic Poisson brackets 
(1) has form 



(la) 


and Jacobi identity is satisfied 



( 2 ) 


□ 

If / < 4, J < 4 and K < 4 then execution of (2) is obvious. 

Let’s consider completely opposite case when / > 4, J > 4 and K > 4. 


We have 
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(2£y r /i r )’ *^Z+4,fcH-4 “1“ (2^jfcr/^r) 5 *^Z+4,i+4 “I - (2 ^kirf^r)' t ^Z+4,j+4 

(2£^>5 r ) (2i£ij iS fJ's) 4~ (2<Sjfc r 5 r ) (26‘^ s /i s ) 4~ (2£^ r 5 r ) (2£ys/x s ) 

4(^ijr^-r/cs 4“ ^jkr^ris 4“ ^kir^rjs)f^s 
4((5^^'s ^is^jk) 4“ (fiji^ks fijs^ki) 4“ ( ^kj^is ^ks^ij')')l^s 0 

If I < 4, J < 4 and AT > 4 then 

Jap' JL,k-\-A 4“ J/3,k-\-4’ JL,a 4“ Jk-\-4,a' Jl,P 
JL,a 4“ Jl,(3 «/fe+4,/5’ JL,a 4“ *^+4,a* Jl,P 

Jk-\-4,p'^ J“y,a 4“ Jk-\-4,a'^ J^,p 0 
If now we have I < 4, J > 4 and jRT > 4 then 


(3) ^aj+4 e /l/,fc+4 4“ *^7+4,A:+4 JL,a 4“ «^A:+4,a JlJ -\-4 

= —Jj+4 

,a' f J"/,k+4 + ■/j+4,fe+4’ Z+4 <^+4,a + <^fc+4 ,a’ ^7^+4 

Jj-\-4:,a 5 ^ t /fc+4,7 «^A;+4,a ,/ ^*^7+4,7 4“ ^7+4,/c+4’ *^Z+4,Q! 

And also let now <a = 0 then expression (3) takes the form 

*/j+4,0 , ^*^fe+4,7 Ac+ 4,0 ^^’+4,7 4“ ^’+4,£+4’ *^i+4,0 

^7+4,0 e /fc+4,s t /fcH-4,0 5 Jj-\-A,s 4“ 2 SiklQl 
^j\^ksrQr Qo^ks^) ^kK^jsrQr Qo^js^) 4“ 2 &iklQl 

kjrQr Qo^kj^) jkrQr Qo^jk) 4“ 2 &iklQl 

2 £jkrQr 4“ 2 SiklQl 0 

Else if a = t then expression (3) takes the form 

Jj-\-A,t 5 ^ e /fc+4,7 Jk-\-4jf^ Jj +4,7 4“ Jj+4, /c+4’ *^i+4,t 

Jj-\-A,t Jk-\- 4,0 *^feH-4,t *^7'+4,0 4” 2<S jklJl-\-4,t 4“ Jj-\-A,t e /fc+4,s '/fc+4,d ^/+4,s 

fijtQk 4“ &ktQj 4“ 2 £jkl\£ltrQr Qo&lt) 4“ £jts\£ksrQr Qo^ks^) ^kts^jsrQr Qo^js') 

3jtQk 4“ ^ktQj 4“ 2<S jkl^ltrQr 4“ ^jts^ksrQr ^kts^jsrQr 

— 2 £jktQ0 ~ e jtkQ0 + ZktjQ 0 

fijtQk 4“ $ktQj 4“ 2 Sjkl^ltrQr 4“ ^jts^k srQr £kts£jsrQr 

fijtQk 4“ &ktQj 4“ 2(fijtfikr 3jr^kt)Qr (&jk&tr &jr&tk)Q_r 4“ (&kj&tr ^kr&tj ) Qr 

fijtQk 4~ &ktQj 4~ 2 SjtQk 2 SktQj &jkQt 4~ $tkQj 4~ &kjQt &tjQk 0 
Lists all possible combinations of indices. 
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8.2 Poisson map from S 3 x M 3 to S0( 3) x so(3)* 

Let’s suppose now that on S 3 x M 3 we have next Poisson structure 


'{Q^,Qu} = 0, fi,u = 0,1,2,3; 

/-j\ ^ {/L> Qo} Qii 

{l*i i Qj } £ijkQk Qo^ij i 

kl/hj/b/ ‘^Eijlfth 

Let’s put in dynamic variables 

(2) 7b = 

Proposition 5. If Poisson brackets in quaternion variables (1) is valid 
then for matrices in form 

(3) Qik (2g 0 T 2q^q k 2,qoQjEjik 

the next Poisson brackets will fulfilled, 

(4) {QijjQkl} 0) {tTj, Qjk} E-ijlQlki {Tq - 77j | EijlT^l 

i.e. Poisson structure on T*SO( 3) is induced from Poisson structure (1) on 
T*S 3 put it in other way T x ^id R 3 is a Poisson mapping. 

□ 

From (1) follow next Poisson brackets 

{ {7h,<?o} = <m*; 

{ t Ti,qjq k } = \£iknqjqn + \£ijnqkq n ~ \qoqjSik ~ \qoqAj’, 

{ t U, q 0 qi} = -\ql8u + + \E an q Q q n \ 

Then 

{ 7 b, Qjk} = 2{7 n, ql}5 jk + 2{7Tj, q 3 q k } - 2£i jk {n i , q 0 qi} 

1 

^2^iknQjQn 

-2 Eijk (~ql5u 


2 ^ijnqkqn ^qoqj^ik 2 qoqkSij 


l l 

+ 2 ^ + 2 £ilnq ° qn 


— 2qoqiSj k + 2 




^QoQi^jk "b £ iknQjQn “b £ijnQkQn QoQj^ik QoQk^ij 

~^~£ljkQo^il £IjkQiQl £ljk£ilnQoQn 

2 QoQi^jk “b £iknQjQn "b £ijnQkQn QoQj^ik QoQk^ij 

~^~£ljkQo&il £ljkQiQl "b £jlk£ilnQoQn 
^QoQi^jk “h £ iknQjQn ~b £ ijnQkQn QoQj^ik QoQk&ij 

-^-^IjkQo&il £IjkQiQl "b (^ij^kn $ik$jn) QoQn 

^QoQi^jk “b £iknQjQn “b £ijnQkQn QoQj^ik QoQk&ij 

~^~£ljkQ() &il £ IjkQiQl “b qoQkSij QoQj^ik 

Hence we have 


(6) 

{vr*, Q jk } = eijkql^u + ZqoQiSjk - 2 q 0 qjS,, 


“b^ iknQjQn "b £ijnQkQn £ IjkQiQl 

(6a) 

{^ii Qjk} ZljkQo&il “b 2(^0 Qi^jk ^QoQj^i 


~^~£ iknQjQn £ jknQiQn “b £ ijnQkQn 


(66) { 7 ri, Q jk } = SijkQlSu + 2 q 0 qiSjk - 2q 0 qj5 ik 

"b i^ikn^jl £jkn$il) QlQn "b £ijnQkQn 

Since value Sikn^ji — £jkn&ii is antisymmetric with respect to ij then it can 
be expressed by means of dual value. 

(7) 2 £mij (^ikn^jl £jkn$il) &kl&mn &km&nl 

i.e. 

(7 Cl) £ikn$jl £jkn&il £ijm\^kl^mn $km$nl) £ijn$kl ~ £ijk$nl 

Substituting (7a) in (6b) we obtain 

{^ii Qjk} ^IjkQo^il ~b ^QoQi^jk ^QoQj^ik 

"b \£ijn^kl £ijk$nl) QlQn “b £ijnQkQn 
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e ljk%$il "t" ZqoQifijk ^QoQj^ik ^ijkQ. "h 2<S ijnQkQn 
^IjkQo^il “t“ 2 QoQn (&jk$in $jn$ik) ^ij fc (1 9o) 2 ^ijnQkQn 

(2^0 1 A: “1“ 2 QoQn (&jk$in &jn$ik) 2S 'ijnQkQn 

(2^0 "h ‘^QoQn^ijm^nkm ~\~ 2 ^ijnQkQn 

£-ijl ((2^0 l)^Zfc “I - ^QlQk "t" ‘2‘QoQn^nkl) 

= Siji ({2ql - 1)5** + 2 q t q k - 2q 0 q n e n i k ) 

Now then 

( 8 ) {^"25 Qjfc} ((2?0 l)5*fc H - 2 QlQ k 2qoq n S n l k j ^ijlQlk 
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9 Symplectic geometry on T*S 3 ~ S 3 x K 3 


Let’s consider Liouville-form on T*S 3 . 

(1) 0 = Hi{6q)\ 8q = dqq~ l = dqq* 

It follows (d(g _1 ) = —q~ 1 (dq)q~ 1 ) 

d© = dni A (dg)* + Hid(5qY = d/ij A (dg)* + /ijd(dgg -1 )* 

= dfii A (dg)* - Hi(dq A d(g -1 ))* = d/i*(dg)* + Hi(dq A g _1 dgg _1 )* 
Let’s consider 2-nd summand on pair of vectors 

dq A g _1 dgg _1 (w, v ) = dq(u)q~ l dq{v)q _1 — dq(v)q~ l dg(n)g _1 

= dgg -1 A dqq~ 3 (u, v) = 5q A dg(w, n) 

Hence 

(2) d© = d/ij(dg)* + iiid(8q) 1 = d/ij(dg)* + /i*(dg A dg)* 

(2a) fl = —d© = —d/ij A (dg)* — /ij(dg A dg)* 

(26) fl = (dg)* A dyUj — /ij(dg A dg)* 

Consider the expression (dg A dg)* by calculating it on vectors u,v 

(dg A 8q) l (u,v) = (dg(n)dg(n) — 5q{y)8q{u)) 1 

= 2(dg(u) x dg(u))* = £ ik i(8q k (u)8q\v) - 8q k {v)8q\u )) 

= Afc(dg fc A 8q\u,v ) 

i.e. 

(3) (dg A dg)* = Awdg fc A 8q l 
Hence 

(2c) H = dg* A d/ij - /i i£ ik i8q k A dg* 
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9.1 Symplectic structure on S 3 x M 3 that is induced 
from SO( 3) x so(3)* 


(1) ^L 5 ° (3) = -d0 T * so(3) = dQ* A dTTi - 7Ti[ 5Q, 5QY = 

Let’s express differentials 5Q by means of dg as well as using dit = |d/i 
We have 

~£ijkSQjk A d-Ki = 5Q l A d-Ki = 2dg* A ^d/p = dg* A d/x* 

and also 

^ ijk^Qjs A 8Q sk (jll : V ) E-ijk^dQjs^Zl'jSQsi c(u) dt^ js{v^)dQ s/c(^)) 

= = £ijk(5Q(u) X d(2(u)) jfc = -2 (SQ(u) x dQ(u))* 

= -8(dg(xx) x dg(w))* = -8 e ijk 5q ] (u)6q k {v) 

= -Ae ijk (5q J (u)8q k (v) - 8q>(v)8q k (u)) = -Ae ijk {8q ] A 8q k )(u,v) 

i.e. 

(2) £ijk8Qjs A SQsk(u, v ) = A dg fc )(xx, u) 

or 

(2a) EijkdQjs A 8Q sk 4 £ij k (5q^ A dg ) 

Thus 

f^can -^ZijkSQjk A d7Tj T ~^i£ijk^Q js A 8Q sk 

= dg* A d/i, - 2n i £ i j k (8q :i A dg fc ) = dg* A d/x* - /x i &y A; (dg- ? A dg fc ) 

that coincides with (2c) of previous subsection. 

Then through direct calculation we proved in representation of right triv- 
ialization next proposition. 


Proposition 6. Symplectic structure on T*!? 3 is induced from symplectic 
structure on T*SO( 3), i.e. 


(3) 


(r 


x 


;id R >)'ei; so(3 > = e£,f; 


(r x iid R .)'a£f° l3) = n™’; 

and this mapping is revealed local isomorphism. 
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9.2 Left- and right-invariant field on T*S 3 ~ S 3 x R 3 

First of all, left and right action on T*S 3 ~ S 3 x R 3 takes the form 


( 1 ) 


Then 


R P (q, v) = {qp, /*); 

Lp(q, p) = (P9,2W _1 ) = (pg,r(p)[/j]); 


i R (q,v) = (q£, o); 

i L (q, I*) = (&, - m£) = (£?, 2 £ x /j); 


From (2) we have next proposition. 


Proposition 7. The vector fields 


(3) 


e* 


A 

’ dni 


form the global basis on T*S 3 ~ S 3 x M 3 . 
Then an arbitrary vector 


(4) 


X e T(S 3 x M 3 ) 


v k £k + Vk 


d 

dpk 


9.3 Poisson structure of symplectic geometry on T*S 3 
We have 


( 1 ) 


Then 


L^ r Q — 0; 
L^l 0 = 0 ; 


(2) = —i^clQ = —+ do i^(~) = d(0(£)) 

and 


ieRto = d(G(q£,o)) = d((ii(q£q l ) i ) = d(v,q£q x ); 
= d(Q(£q, 2£ x p)) = d(fuC) = d{n, £); 
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2-nd line in (3) means that Hamiltonian field corresponding to dynamic 
variable (/x, £) is £ L . 

Then 

(4) {(/l€>, = 2(?7 = 2(/w,£ x ry) 

(4a) {/x, (pt,r 7 )} = 2r/x pt 

(46) {/ii,p fc } = 2e ifcZ /i/ 

Analogically 

(5) {^(,x,0} = %g Q = (£<?)“ 


or 


(5a) {<?, (/x, 0} = & = ~(£, <l) e o + <?°£ + £ x q 

Hence we have 


m 

\{<f, (/*,£» = + (£ x <?)b 

(5c) 

- (£ X q)* 

(5 d) 

f{/L,<? 0 } = g»; 

= -g°<5f + e ik iq l ] 


Now let’s consider Poisson brackets for functions that depend only on 
variables q but not /i. 

For basis quaternions from (3) we obtain 

(6) AlQ = dfik 

Let coefficients v k are arbitrary functions (q, /x). Then we have 

(6a) i^ v k^Q = v k dfik 
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Let 


(7) 


X = v k ~e L k + Vk ^- e T(S' 3 x M 3 ) 

o/ifc 


arbitrary vector. 

From (6a) and expression 


( 8 ) 

we have 


Q = (. Sq) l A dni - HiE ik M k A Sq l 


(9) i x &. = v k d/d k + q k 5q k 

In equation for Hamiltonian field Xp corresponding to Hamiltonian F(q) 

(10) ix F M = dF 

as it is follows from (9) the held Xp can not contain components related to 
basis vectors e k i.e. 

(11) X F = q k -— 

Olik 

If F and G are two such functions then from 


(12) {F,G}(«) = n(^),fc(a;)) 

it should be 


(13) {F{q),G(q)} = a 

and hnally we have next proposition. 

Proposition 8. The Poisson structure of symplectic geometry on T*S 3 
is dehned by next Poisson brackets 

({F(q) 1 G(q)} = 0; 

(1 4 ) \ {?, (Mil)} = £9; 

(Ab 7 ?)} = (^1 t^ 7 ?]) = 2 (/b£ X T]); 

Remark 1. Poisson brackets (14) are revealed coordinateless form of 
Poisson brackets (1) §8. 
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Remark 2. So that get Hamiltonian field for variables qo,qi from 


(15) 


>x.P- = Vk&q k 


one can use the matrix of right action in the space of quaternions 


( 16 ) 


Rb — 


b° 
b 1 
b 2 
b 3 


-b 1 

b° 

-b 3 

b 2 


-b 2 
b 3 
b° 
— 6 1 


—b 3 

-b 2 

b 1 

b° 


Then we get 


(17) r/ 0) = 

1 

1 1 

to M 

, 'n {1) = 

i 

CO o 

_1 

, ?7 (2) = 

— q 3 

qO 

, t 7 {3) — 


— q 3 


1 

CN 

1 


1 
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10 Hamiltonian equations of motion for a 
rigid body in quaternion variables 


10.1 Poisson structure of a rigid body dynamics in 
quaternion variables in inertial frame of reference 

As shown in [11] right trivialization of cotangent bundle of T*SO( 3) and T*S 3 
for description of a rigid body dynamics has the following advantage: as far as 
translational and rotational degrees of freedom are separated in inertial frame 
of reference; then Poisson structure can be represented as direct product of 
structures that are refer to these degrees of freedom of rigid body. 

Therefore, adding Poisson structure of subsection 8.1 by Poisson brackets 
for translational degrees of freedom we obtain 

{xi,x k } = 0 , {x i ,q fl } = 0 , {xi,Hj} = 0 ; 

fe.,^} = 0 > {Pi^j} = °; 

{xi,Pk} = $ik, i,k = 0,1,2; 

< = 0, fi,u = 0,1, 2, 3; 

{ h i j Qo } Qii 
{Pit Qj} £ijkQk 

= 2 Eijlfll, 

Poisson structure in quaternion variables in mixed 
frame of reference 

If we have Poisson brackets of previous subsection and system Hamiltonian 
we can write Hamiltonian equations of motion. 

But, as well known, expression of kinetic energy of a rigid body is much 
more simpler in the body frame rather than in inertial frame of reference. 

Therefore, it makes sense to conduct the canonical transformation of gen¬ 
erators of our Poisson structure, namely intrinsic angular momentums of the 
body 7 — so that they represented to the body frame. 

n = Q-'M; 

P - p, X — x, Q-Q 



( 1 ) 


10.2 
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In quaternion variables 1-st line will be as follows 
(2) M = q- 1 fiq = Ad*[fi\ = 


or 

(2a) Mj = QrjUr 

That was revealed previously (see (8) in subsection 8.2) 


(3) 


{A Hi Qjk} ^E-ijlQlk 


Therefore 


{ALjMj} }j^ii Qrjf-^r} QrjifJ'ii fJ'rJ T Qrj}f^r 




0 


i.e. 

(4) {^t,Mj-} = 0 


Consider also 


{Mj, Qjk / {Qril^ri Qjk} 


Qri{l^n Qjk} rktQriQtj 


Let’s use identity 


(5) 

from it implies 
(5a) 


^rstQriQtjQ sk £ikj 


^-rstQri Qtj Qkl^ilj 
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Hence 


Qjkf 2 SrktQriQtj ^Qkl^ilj 2 SijlQkl 


( 6 ) 


{Mj, Qjkj ^E-ijlQkl 


Now consider 


i.e. 


{ M, Mj} — {M j, QrjUr } — {Mj, Q r j } fi r + Qrj'{Mj, /i r } 
{Mj, Qrj}^r ^^ijlQrl ^IJ /M/ 


( 7 ) {M i ,M J -} = - 2 e iiI M, 

Let’s obtain Poisson brackets between moments and quaternion dynamic 
variables 

{Mj, go} {Qril^ri go} Qri{^ri go} 

QriQr (( 2^0 "P 2 Q r Qi 2 QoQl^-lri) Qr 

= ( 2 go - !)& + 2 Q 2( li = Qi 
i.e. 

(8) {Mj, g 0 } = qi 
Now consider 


{ 1 VI ii Qj} {Qri/J'riQj} Qri{H j nQj} 

((2# 0 d" 2 (JrQi 2q , q Ql^lri) rjsQs Qo^rj ) 

((2^0 d" ^QrQi ^QoQl^lri) £rjsQs 

Qo ((2% d" ^QrQi 2(/o Ql^lri) &rj 

(2^0 1 'j^ijsQs ^QoQlQs^lri^rjs 

-<7o(2<7o - l)£y - 2<Mj% + Zqlqi£iji 
= {2ql - l)£ijiqi - 2qlq i £ ijl 
-?o(2go “ l)$y - 2<Mj& 
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~\-‘^QoQlQs(^jl&is ) 

= Sijiqi - <?o(2<?o - l)^ii - 2g 0 ftft + 2g 0 g i g i - 2g 0 <7% 

= -£yi9z - <?o(2<7o - !)% - Zqoq 2 $ij 
= Sijiqi ~ <?o(2<?o + 2<Z 2 - l)5jj = Sijiqi - q 0 $ij 

i.e. 

(9) {Mj, qj} = —qo$ij - Eijiqi 

Hence finally have next proposition. 

Proposition 8. The following Poisson structure in quaternion variables 
for rigid body dynamics in mixed frame of reference is valid. 


( 10 ) 


'{fti,ft/} = 0, n,v = 0,1,2,3; 

{Mi, go} = ft; 

{Mi, qj} = ~q 0 S i:j - Eijiqi ; 

, {Mi, Mj} = -2ei j7 M,; 


10.3 Motion equations of a rigid body in mixed frame 
of reference 

Let’s consider the Hamiltonian of rather general form for Hamiltonian dy¬ 
namics of a rigid body in quaternion variables. 

(1) H((x,p),(q, M)) = p 2 + T spin ((x,p), (q, M)) + V(x, q) 
where 

1 1 /M 2 M 2 M 2 \ 

(2) T spin ((*, p), ( q , M)) = + -j4 j 

In mixed frame of reference arithmetic vectors (x,p) (i.e. translational 
degrees of freedom) are components of the corresponding physical vectors in 
inertial frame of reference. And arithmetic vector M (rotational degree of 
freedom) are components of the corresponding physical vector in body frame; 
I — constant tensor of inertia in the body frame. 
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Using base Poisson brackets in mixed frame of reference 


(3) 


{xi, x k } = 0, {Xi, q^} = 0, {xi, M j} = 0; 
{P%,%} = 0) {Pii Mj} = 0; 

{xi,p k } = S ik , i,k = 0,1,2; 

{Qn,Qv} = 0, P," = 0,1,2,3; 

{M^ g 0 } = Qu 

{Mi, qj) = -q 0 Sij - Eijiqi ; 

{Mi, Mj} = -2^ 7 Mi; 


for the Hamiltonian (1) we obtain the next Poisson brackets 


Stand for 


(4) 


±i = {xi, H} = Xi, -^—p 2 1 = — pi 
2m m 


pi = {p^ H} = {pi, V(x, q)} = 


dV 

dxi 


Q i = H = 2-^- = -(I-^ikMk 
dTTi dM ,■ 2 V ’ 


(4a) 


1 2 /i 2/ 2 2/3 

dH 

go = {g 0 , H} = — {go, Mi} 

= {g 0 , Mi} = -^Oigi = ~{Q,q) 


i.e. 


( 6 ) 

in what follows 

1 


go = {go,#} = --(O ,q) 


dH 

Qi = {ft, H} = (gi, M fc } 


{ft, Mfc} ^(<?o4; £kuqi) ^qo^i 4~ ^^kuQi^k 
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— 2 < ?0^i ^ikl^kQl 


i.e. 


(5) 


q,, = Ui, H} = ig 0 hh - ^Eiki^kqi 


or 


(6a) 


q = ^q o n - ifi X q 


in what follows 

r)M r)H r)M 

Mi = {Mj, H} = —— {M,, g 0 } + w— {Mi, g*,} + „ {M* 

ago oqk aM^ 

dV dV 1 

TJ q% ~h TJ ( qo^ik Sikiqi) ~t~ ( 2£ifc/M/) 

<9g 0 dq k 2 

dV dV dV „ ,, 

77 g* g0 77 E-ikll^ qi *fcM; 

dg 0 oqi dq k 


i.e. 


(7) 


• dV dV dV „ ,, 

M* = —gi - g 0 ^- £iki-^—qi ~ eikiiliMi 

dq 0 oqi dq k 


or 


(7a) 


M = ^-q - q 0 d q V - d q V x q - flxM 
dq 0 


Now then 


( 8 ) 


X =™P> 

P = ~d x V] 

< g 0 = {g 0 , H) = -|(0,q); 
q = |g o n - |0 X q; 

= -n X M + gq - qodqV - d q V X q; 


Let’s lump together 4-tli and 5-tli line in (8). 


(9) 


1 . 1 ^ 1 
Q = ~ g q)e 0 + -g 0 n - -S2 x q 



Mfc} 


42 





(9a) 



Consider the quaternion 

(10) V (9 V = e 0 d qo V + d q V 


Then 

q-'V^V = q'V®V = (g 0 e 0 - q ) ( e 0 d qo V + d q V) 

= (q, V^V)e 0 + q 0 d q V - d qo Vq - q X d q V 

i.e. 

d qo Vq - q 0 d q V + q x d q V = (g, V (g) C)e 0 - <tV 9 V 
Therefore 

(ii) M = ^nxM-(f 1 v (,) H?,v (, V)) 


( 12 ) 


P = -d x V; 

i -i i 

q = 5 gO —> q i q= 

M = (g^VWy - (g, V (9) H)) ; 


Proposition 9. The motion equations of a rigid body dynamics for 
Hamiltonian (1) in quaternion variables has the following algebraic form 


(12a) 


P = ~d x V] 

q = |gO —> g _1 g = |0; 

M = -n x M - 9= (g-'V^H) ; 


where 5s(a) — imaginary component of a quaternion a. 
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